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Abstract. 

Photonic bandgap in holographic grating manifests itself as phase-sensitive 
birefringence. Phenomenological theory of the effect is presented. 



PACS numbers: 42.25.Lc, 42.40.Eq 
1. Introduction 

Microscopic mechanisms behind the formation of holographic grating are complicated. 
Respective nonlinear system of differential equations for the kinetics of the process, 
which includes many parameters that are used to characterize the photorefractive 
material, can be solved through an intricate numerical calculation only, along with 
making a series of simplifying assumptions. Therefore, it is interesting to construct 
simpler phenomenological models to cope with the problem. 

A simpler approach can be developed towards the non-equilibrium steady states 
(NESS) that are reached when the diffraction efficiency of the grating changes no more 
under stable irradiation [lj. This is a NESS that will be characterized here by a single 
mode of the electromagnetic field (the writing mode), and the respective modulation of 
the electric permittivity (the photorefractive response). The two latter quantities are 
interdependent, and should be obtained in a self-consistent manner in the theory. (One 
may think of some functional, which has extrema at NESS "points" of the system driven 
by irradiation.) 

From the viewpoint of nonequilibrium thermodynamics, the driven open system 
that is represented here by the electromagnetic field coupled to the photorefractive 
medium, slowly evolves towards a NESS under fast external driving. The two time- 
scales differ in a dozen of orders in magnitude, and it is a natural phenomenological 
assumption that the modulation of the permittivity of the initially homogeneous medium 
is determined by the NESS mode itself. Based on this assumption is the discussion 
below. 

An interesting feature of locally isotropic photorefractive medium - phase-sensitive 
birefringence effect, that is a manifest of the photonic bandgap and Bragg reflection, 
shows up. 
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2. Exactly solvable model 

Vector nature of light is at the heart of the photorefractive effect. Fortunately, the two 
polarization modes can be treated separately, in the properly symmetric arrangements. 
In view of the holographic grating recording, of major interest are the modes of the 
interference pattern. Note that their observation can be perfectly "tuned in" due to 
Bragg reflection. By assuming that the active NESS mode is the electric component of 
the so-called E-wave, 



Maxwell's equations are reduced to the second-order differential equation for the mode 
profile, F(z): 



for details see [2]. F(z) is referred to below as the "NESS mode". It is assumed that the 
in-medium mode inherits i) the antiperiodicity of the driving field, F(z+X z /2) = —F(z), 
and, similarly, ii) the reflection symmetry for the mode profile, F(—z) = F(z), or 
F(— z) = —F(z), and, consequently, the eq. © has to be solved under antiperiodic 
boundary condition for the even or the odd mode. 

Perhaps the simplest model assumption for the NESS permittivity is that that its 
modulation follows the intensity of the NESS mode directly; then the transparent and 
nonmagnetic medium of the model, with periodic inhomogeneity in the z-direction, is 
described by its permittivity at a given frequency k of the mode by 



Fq is the amplitude of the NESS mode, and it is assumed that the t]Fq is relatively small 
number: the modulation —i]FQ<e(z)—e<r]FQ , typically, is weak in holographic gratings. 
An expression similar to (j3J) is used in nonlinear electrodynamics to study self-defocusing 
of the electromagnetic waves (see [2]). Here its meaning is different: we consider the 
linear electrodynamics problem, and the "nonlinear" term describes the mode profile 
in the NESS. In view of the localization properties of electrons in photorefractive 
materials, which are known to be good insulators, the local approximation in (j3J) seems 
to be reasonable: it is a natural phenomenological assumption that the modulation of 
the permittivity of the initially homogeneous medium, in the first approximation, is a 
quadratic function of the NESS mode amplitude. 

After substituting e(z) as a function of the "brightness" of the mode, e(z) — > 
e — 2i]F 2 + i]Fq, the linear equation (j2J) takes its nonlinear form: 



E y (z, x,t) — F (z) sin{k x x — tut), 



(1) 



F" + {kle{z) - k 2 x )F = 0; 



(2) 



e{z) = e-2r l F(z) 2 + r ] F 2 ; 



(3) 
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Figure 1. Permittivity modulations for < m < 0.02. Plotted is the right-hand side 
of the eq. ©. 



see textbooks on classical mechanics. The mode of interest is, naturally, that of winding 
number 1. Resulting dispersion relation has the form: 

\l(kle-kl) = (AK(m)f = 
(2vr) 2 '- m 



m 



32 



+ 0(m 3 )), 



(4) 



in contrast with the relation XHk^e — k 2 ) = (27r) 2 for an isotropic and uniform medium. 
(We have avoided using the notation " k z " along with the " k x " to stress the fact that the 
former is actually a quasimomentum.) K(m) is the complete elliptic integral of the first 
kind, and its modulus, m, is related to the parameters of the problem by the following 
equation: 



m 



V*3 



k 2 



h 2 i= — h 2 ' 

n e n x 



The NESS mode has the form 

AK (m)z 



F(z) = F sn 



A, 



with the Jacobi elliptic function sn. The respective profiles of the permittivity 
modulation normalized to the effective NESS mode, 

\l(k*e(z) - kj) (AK(m)) 2 ^ 



(2 



7T 



2m I sn 



(2tt) 2 
4K(m)z 



rn 



+ m 



(5) 



A 2 

are shown in the figure [U 

The eq. (|2J) with e(z) defined using (|3j) reduces to the well known equation of Lame. 
With the brightest point of the interference pattern at z = the Jacobi sn in (jlj) and 
(J5J) will be replaced by the Jacobi cd. For references on Jacobi elliptic functions and 
Lame equation see [3]. 
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Remarkably, the inverted brightness (F*(z)) 2 = F 2 — (F(z)) 2 mode also exists as 
an eigenmode for the eq. <^j, 

AK(m)z 



F*{z) = F cn 



in 



the cn will be replaced by the y/1 — msd when the brightest point of the interference 
pattern is at z = 0. The dispersion relation for these latter modes has the form 

X 2 z (k 2 e-kf) = (l-m)(AK(m)f = 

(2vr) 2 (l-^-^ + 0(m 3 )), (6) 
and, therefore, the difference 

2 

777 

(2n) 2 (m+ — + 0(m 3 )) (7) 

is small as expected. The picture is sensitive to the relative phase, and the dispersion 
relations (jl]) and (J6]) can be contrasted with those for the conventional birefringence 

case, 

>?Me - kj) = (2tt) 2 , \ 2 z (k 2 e e - kj) = (2vr) 2 ^, 

with the indices o and e denoting the ordinary and extraordinary waves, respectively. 

The two solutions found belong to the band edges of the spectrum of the 
Hamiltonian associated with the periodic modulation of the e(z); the bandgap here 
manifests itself as a phenomenon resembling birefringence: the above expression (El) 
describes phase- dependent birefringence; joint dispersion relation in parametric form 
is given by (j4|) and (jSJ), and the small parameter m is accounted for by (J7j). With the 
latter equation, experimental observation is straightforward, by measuring the difference 
between k* and k x . Regarding such an experiment, as well as the symmetries of the 
problem, the remarks in the following section seem relevant. 



3. NESS, two basic geometries, and translation invariance 

As regards the (relative phase-sensitive) "Bragg birefringence" effect mentioned at the 
end of the previous section, two setups shown in the figure [2] have to be considered. 
Full translation invariance in the bulk of the photorefractive medium is conserved along 
the horizontal axis only (the x-axis). With the boundaries of the plate parallel to the 
vertical axis (the z-axis), externally controlled is the quantity A z , and the quantity k x 
- in the next setup, with the boundaries parallel to the x-axis. Also notice that the 
parameter e that enters the joint dispersion relation flj]) and ([6]) is " renormalized" by 
the NESS grating, and actually also depends on the m, as a self-consistent quantity. 
Therefore, to use the dispersion relations, the small parameter m has to be determined 
first, by using the relation ((?]). This can be done, in this example, due to the reflection 
symmetry: by externally exciting the two modes, F and F*, separately. 
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Figure 2. Two complementary setups for holographic grating recording. The two 
incident plane waves dictate either X z (on the left) or k x (on the right). 

Interference pattern created by two crossing plane waves in a homogeneous 
transparent medium - the intensity of the superposition mode of the electromagnetic 
field, is considered to be the "driving" field. For instance, the superposition mode (the 
two incident waves have opposite phases at the origin of the coordinate system) 

E y sm(k z z) sm(k x x — tut) 

creates a one-dimensional interference pattern with intensity profile \{E y sm(k z z)) 2 , and, 
similarly, the intensity profile is \{E y cos(k z z)) 2 for the mode 

E y cos{k z z) cos(k x x — ut) 

(for coinciding at the origin phases). As regards the interpretation of the ([3]) in view 
of nonequilibrium processes for the E-wave, the entropy production in transparent 
materials is actually the Joule heat, and the modified E 2 is at place in the equation. 

The period of the interference pattern is X z /2, and the bright and the dark strips 
alternate at half that length. Translation by A z /4, 

F(z)^F(z)=F(z + X z /4), (8) 

results in swapping them (recall that A z /2-antiperiodic modes are considered, and, 
therefore, the direction of the shift could be reversed as well). For the brightness of 
the interference pattern in a homogeneous medium, the following symmetry property 
holds: 

F(z) 2 = F 2 - F(z) 2 . (9) 

For the NESS, however, full translation invariance of the medium is broken due to the 
electric permittivity modulation: the "polaritonic" nature of the NESS mode should 
not be symmetric neither with respect to the bright and the dark, nor with respect to 
the A 2 /4-shifts. This is clear with the quartet of the modes obtained in the previous 
section, which solutions can be denoted as F, F, F*, and F* in the view of the ([8]). Both 
pairs of the modes, F, F*, and F, F*, are described by a pair of the wavenumbers k x , k*, 
and the symmetry ([9]) is broken. 

The symmetry with respect to shifts by A 2 /4 (see the definition (jSj)), obviously, 
is not for the homogeneous medium uniquely. With full translation invariance in the 
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Figure 3. Symmetric modulations: < m < 0.1. 

z-direction broken, there can exist modulations that preserve that symmetry. An entire 
family of such is presented by a symmetric modification of the eq. (J3]): 

\ 2 z (k$e(z) - k 2 ) (AK(m)) 2 ^ 



(2tt) 2 (2tt; 



2 



m 3 , n .2 F 2 + F 2 



2 4 v \ \ F, 



2 



The respective profiles of the permittivity modulation are clear from the figure |3j 
plotted is the right-hand side of the above equation. The problem is equivalent to the 
so-called associated Lame potentials [I]. In this symmetric case, the two modes have 
the same k x : no Bragg birefringence is present. The modes, explicitly, are: 

F (l — m) 1//4 sn (u\ m) ~ F cn(u\m) 



dn(u\m) ^dn(u\ m) 



with u = iK ^ z . The Wronskian W(F, F) (defined as FF' — FF') is constant, because 
of the "degeneracy" mentioned above. This time the symmetry with respect to the 
inversion of the brightness does not hold, in contrast with the model presented by the 
eq. ([3]): the respective mode F*(z) such that (F*(z)) 2 = Fq — (F(z)) 2 does not exist as 
an eigenmode for the given modulation of the permittivity; this is quite natural in the 
view of the degeneracy present in the second-order equation. 

Gratings of this type can be of interest for the optically nonlinear photorefractive 
materials in the setups appropriate for the higher harmonic generation. 



4. Concluding remarks 

In the view of the variety of possible mechanisms of the holographic grating formation 
in photorefractive media, it should be noted that even in a single material there may 
exist different mechanisms, due to the variations in the irradiation techniques (e.g., with 
pulsed irradiation, as compared to the cw irradiation, by varying the pulse duration and 
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the intervals in between vs electronic relaxation times). It means that engineering the 
models is also possible. 

In the case the holographic grating is recorded in a plate with /^-controlled setup 
shown in the figure [2j its thickness I comes into play: because of reflections, it must 
match the phase of the interference pattern such that 21 = n\ z . Similarly, in the case 
of using a mirror for recording, such considerations could help for better tuning the 
process, by slight variations of the angle of incidence and, if possible, the frequency of 
the incident beam. 

In conclusion, it was discussed how the permittivity locally-isotropic modulation- 
induced anisotropy appears, leading to the Bragg birefringence effect, in the process 
of the recording of the holographic grating. The effect can be of interest for quantum 
devices, due to its phase-sensitivity. 
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